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The forces between charged macromolecules, usuaUy given in terms of osmotic pressure, are highly 
affected by the intervening ionic solution. While in most theoretical studies the solution is treated 
as a homogeneous structureless dielectric medium, recent experimental studies concluded that, for a 
bathing solution composed of two solvents (binary mixture), the osmotic pressure between charged 
macromolecules is affected by the binary solvent composition. By adding local solvent composition 
terms to the free energy, we obtain a general expression for the osmotic pressure, in planar geometry 
and within the mean-field framework. The added effect is due to the permeability inhomogeneity and 
nonelectrostatic short-range interactions between the ions and solvents (preferential solvation) . This 
effect is mostly pronounced at small distances and leads to a reduction in the osmotic pressure for 
macromolecular separations of the order 1-2 nm. Furthermore, it leads to a depletion of one of the 
two solvents from the charged macromolecules (modeled as planar interfaces). Lastly, by comparing 
the theoretical results with experimental ones, an explanation based on preferential solvation is 
offered for recent experiments on the osmotic pressure of DNA solutions. 



I. INTRODUCTION 

The interactions between charged macromolecules im- 
mersed in aqueous solutions are of great importance 
in biology and material science. Because of their rel- 
evance to colloidal suspensions and biological macro- 
molecules, the forces between charged objects mediated 
by electrolytes have been the focus of numerous stud- 
ies [H, 0, i, i, H, H, 0, B i- Within the context of 
the so-called primitive model, the solvent is modeled 
as a homogeneous dielectric medium [13, El, [13 ^if- 
fecting the system only through the dielectric constant 
that acts to reduce the strength of the electrostatic field. 
More recently, theoretical and experimental approaches, 
studying interactions between charged macromolecules, 
have been extended to also treat binary solvent mixtures 

The thermodynamics of binary solutions is well under- 
stood and has been described in detail in many textbooks 
However, there are still open questions concerning 
the behavior of binary solutions in the presence of other 
degrees of freedom, such as dissolved ions and external 
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electric fields. These additional couplings are relevant to 
a broader spectrum of applications, extending from ma- 
nipulation of microfluids JJi, l2fl, l21j, [22i] to biologically 
motivated problems such as protein stability and con- 
formational changes [2^, [l^, [5^, [2^, For example, 
in recent experiments, the transition of a DNA molecule 
from elongated coil to compact globule was found to de- 
pend on the addition of another polarizable solvent to 
the aqueous solution , suggesting that the interaction 
between DNA segments is modified by the presence of 
this additional solvent. 

The effects of adding cosolvents to aqueous solutions 
of charged macromolecules can be quite pronounced. In 
fact, one of the common ways to precipitate DNA in- 
volves adding an excess amount of ethanol to the aque- 
ous solution, which counteracts the repulsion between 
charged DNA strands [3, [1, @1- This effect has been 
commonly attributed to the change in solution dielectric 
constant. However, studies over the last decade convinc- 
ingly demonstrated that alcohol changes the disjoining 
{i.e., the interaction) pressure between DNA strands to 
a much greater extent than would be expected from the 
direct change in the dielectric constant. This added effect 
that goes beyond changing of the dielectric constant has 
been explained in terms of the preferential exclusion of 
alcohol from the vicinity of interacting DNA strands 0] . 



These studies further demonstrated that because alcohol 
exclusion causes an additional osmotic pressure difference 
between the bulk solution and the concentrated DNA 
phase, DNA strands are pushed even closer together. 

Two distinct features prevail when trying to model 
ions immersed in a binary solvent mixture within the 
standard Poisson-Boltzmann (PB) theory. First, the 
disparity between the solvent permeabilities leads to a 
dielectrophoretic force. The ensuing force acts on the 
solution, attracting the high permeability solvent (e.g., 
water) component toward the charged macromolecular 
surface and, at the same time, depleting the lower per- 
meability one (e.g., alcohol). As a result, the solution 
becomes inhomogeneous and a permeability gradient is 
created in the vicinity of the charged interface, where the 
system favors the higher permeability component that 
can better screen the electrostatic field while excluding 
the low permeability solvent away from charged inter- 
faces. The second important feature is the chemical (non- 
electrostatic) preference of the ions for one of the two 
solvents. The dissolved ions effectively drag with them 
a solvation shell preferentially enriched in one of the sol- 
vents, thus repelling the second. When attracted to the 
oppositely charged surfaces, the dissolved ions thereby 
change the composition of the vicinal solvent. These two 
effects can enhance or compensate each other. In this 
work, we treat only the case where the two effects act 
synergistically to mutually enhance each other. 

Previous theoretical works describing the effects of bi- 
nary solvent mixtures dealt mainly with systems close 
to their critical point. Tsori and Leibler investigated 
the change in the phase transition temperature due to 
dielectric inhomogencity and preferential solvation p^ . 
while Onuki and Kitamura investigated corresponding 
surface tension and the ionic distribution near an inter- 
face d, Ell- To contrast and compare, in the present 
work, we focus on binary solution systems in the sin- 
gle phase region and away from the coexistence region. 
Moreover, contrary to previous works, our main interest 
is the effect of the dielectrophoretic force and preferential 
solvation on the pressure (or forces) between two equally 
charged objects, such as a pair of charged DNA strands. 

We model the system by delimiting ourselves to the 
simple planar geometry for two interacting macromolec- 
ular surfaces. Some experimental setups apply directly 
to this geometry and even for more complex setups, our 
model captures the essential physics of coupling between 
the binary solvents and mobile ions. 

In what follows, we present a model where the ionic 
densities, the solvent relative composition, and the elec- 
trostatic potential are all continuous functions of the lo- 
cal position. We derive a set of coupled differential equa- 
tions relating the various degrees of freedom at thermal 
equilibrium. Furthermore, we derive a general expres- 
sion for the local pressure in the form of a modified con- 
tact theorem and provide proof that it is spatially ho- 
mogeneous. This allows us to reduce the corresponding 
Poisson-Boltzmann equation to a first order differential 
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FIG. 1: Schematic illustration of the model system. The two 
plates residing at z = ±_D/2 are charged with surface charge 
density a < 0. The two solvents are represented by circles 
denoted A and B, with dielectric constants ea and eb < £a, 
respectively. The counterions and the high dielectric solvent 
(sa) are attracted to the plates. 



equation that greatly simplifies the numerical problem. 

Our numerical and analytical results focus on the solu- 
tion mixtures where the low permeability solvent (alco- 
hol) has a small concentration compared to the other sol- 
vent (water). First, we examine the influence of the vari- 
ous parameters on the composition profiles. We find that 
the deviation of the solvent composition profile from its 
average (bulk) value can lead to large deviations from the 
regular Poisson-Boltzmann theory predictions, especially 
regarding interactions between charged macromolecular 
surfaces. We investigate in detail the pressure depen- 
dence on the interplate separation and its sensitivity to 
controllable parameters, such as salt concentration and 
average solvent composition. Finally, we show a com- 
parison between our pressure profiles and the relevant 
experiment on DNA Q. 



II. THE MODEL 

In the model considered here, ions are immersed in a 
binary mixture of two solvents confined between two pla- 
nar charged interfaces. The two surfaces carry homoge- 
neous surface charge densities, cr < (see Figure. [T]). The 
model is formulated on a mean-field level but modifies the 
regular PB theory in two important aspects. First, the 
volume fractions of the two solvents, 0a and (Jjb = ^ — 't'A, 
are allowed to vary spatially. Consequently, the dielec- 
tric permeability of the binary mixture is also a function 
of the spatial coordinates. In the following, we assume 
that the local dielectric response e(r) is a (linear) compo- 
sitionally weighted average of the two permeabilities sa 
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and eb- 



= 0A(r)£A + 0B(r)eB , 



e(r) = £o - </'(r)£r , 



(1) 



(2) 



where we define = ^b, eo = ea and £r = £a ^ £b- 
This Hnear interpolation assumption not only is com- 
monly used but also is supported by experimental evi- 
dence [2^, HO] ■ Note that the incompressibility condition 
satisfies 0a + 0b = 1 , meaning that the ionic volume 
fractions are neglected. 

Apart from long-range electrostatic interactions be- 
tween dissolved ionic species, we also consider the case 
where short-range interactions make an important con- 
tribution to equilibrium properties. Consequently, pair- 
wise short-range interactions between all constituents 
contribute additional terms to the total free energy and 
modify their equilibrium distributions, as will be elabo- 
rated below. 

The equilibrium properties are derived within the 
mean-field framework. Thermodynamic equilibrium is 
obtained by minimizing the (grand canonical) thermo- 
dynamic potential, G = J d^r g{r), leading to a gener- 
alized PB equation that also contains the contribution 
of short-range interactions. The force equilibrium lead- 
ing to interactions between the confining surfaces is then 
obtained from the first integral of the PB equation and 
can be reduced to a surface-normal term of the gener- 
alized stress tensor evaluated at the bounding surfaces. 
The spatial profiles of the two solvents and ions, as well 
as the equilibrium forces, can be obtained from a varia- 
tional principle of the thermodynamic potential. 



A. Free Energy 



We write the bulk free energy as a sum of four terms: 



dV[/,(r) + /i(r)+/„,(r) + /,(r)] 



(3) 



where the free energy density / = /c + /i + /m + fs and V 
is the total volume. The first term is due to electrostatic 
interactions between the ionic species mediated by the 
dielectric medium and characterized by the spatially in- 
homogeneous dielectric function e(r). For simplicity, the 
dissolved ions are assumed to result from a completely 
dissociated (1:1) monovalent salt. In this case, the elec- 
trostatic term, /<,, is given by 

./c = -4^(VV')' + e(n+-n_)^, (4) 

where ip{r) is the electrostatic potential, e is the electron 
charge, and ?T.-|-(r) are the number densities (per unit 
volume) of the monovalent co- and counterions. Note 
that the first term implicitly couples the electric field. 



E = — VV', with the solvent composition, 4>{r), via the 
spatial dependence of the dielectric response, e(r), on 
local composition as was defined in eq [21 This is the 
dielectrophoretic term mentioned previously that favors 
a higher local dielectric constant (lower 0(r)), and causes 
attraction of water to the charged surface. This is one of 
the two sources of the composition inhomogeneity in the 
model. 

The entropy of ion mixing constitutes the second term, 
/i, given by 



/i = fceT 



(log(n+a^) - 1) + n- (log(n_a^) - l) 



(5) 

where fcs is the Boltzmann constant, T is the tempera- 
ture, and is the molecular volume. The third term, 
/m, accounts for the binary- mixture free energy given in 
our model by regular solution theory: 



/m — 



[0 log + (1 - 0) log(l - 0) + x0(l - 



(6) 

The first two terms represent the solvent entropy of mix- 
ing, while the third represents the bilinear short-range 
interactions between the two solvents. The interaction 
parameter, Xj is dimensionless (rescaled by fesT). Note 
that we took the same molecular volume for both A 
and B components. In general, this is not a serious defi- 
ciency and can be easily amended, if necessary. 

The fourth term, fs, in the free energy originates from 
the preferential interaction of the ions with one of the 
two solvents as described in the Introduction. We assume 
here that this preference can be described by a bilinear 
coupling between the two ion densities n± and the rela- 
tive solvent composition 0, which is the lowest order term 
that accounts for these interactions (higher order corre- 
lation terms can be added in a systematic way). The 
preferential solvation energy, fs, is then given by 



fs = kBT {a+n+ + Q!_??._) 



(7) 



where the dimensionless parameters a± describe the sol- 
vation preference of the ions, defined as the difference 
between the solute (free) energies dissolved in the A and 
B solvents. The free energy fs corresponds closely to 
the Gibbs energy of transfer from one solvent to another, 
as is described in detail in refs [3l|, Is^]. This bilinear 
coupling represents a second source of composition inho- 
mogeneity. Namely, a density profile of the ions n±(r) 
(a diffusive layer near a charged object) forces a corre- 
sponding solvent profile, 0(r). 

To all these bulk terms, one must add a surface term, 
describing the electrostatic interactions between charged 
solutes and confining charged interfaces. This surface 
term is given by 



(8) 



where ■0s is the electrostatic potential evaluated at the 
bounding surfaces and depends on the surface charge 
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density a (charged groups per unit area) and surface area 
A. Note that the charged surface is described by a uni- 
form charge density, a. In a more refined model, nonelec- 
trostatic interactions, such as preferential adsorption of 
the two solvents, and finite ion effects could be included 
as well. 

The total free energy is then written as a sum of the 
bulk and surface terms 



d r/(-0, n±, (/)) + d) d r ea^ps 
V J A 



(9) 



In the grand-canonical ensemble the corresponding ther- 
modynamic potential is given by 



/i+?i+(r) + /i_n_(r) + /i^- 



(r) 



g(r) = /(r) - fceT 

where n± and /i^ are the dimensionless chemical poten- 
tials coupled to the ionic densities n± and the relative 
solvent composition 0, respectively. 

In thermodynamic equilibrium, the spatial profile of 
the various degrees of freedom characterizing the system 
is obtained by deriving the appropriate Euler-Lagrange 
(EL) equations via a variation principle of the thermo- 
dynamic potential, eg 1101 The EL equations are then re- 
duced to four coupled differential equations for the four 
degrees of freedom, V'(r), n±{r) and (f){r): 



V • — VV' + e(n I - ?i_) = 
\47r ' 



± 



+ log(n±a^) -I- Q!±0 - ^± = 



(11) 



(12) 



log 



;(V^)2+X(1 
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-I- (a+n+ + a_n_) — /i^ = . 



(13) 



At the charged interfaces, an additional equation stems 
from the surface term of eq[H] 



= 



n • VVj 



(14) 



where n is the unit vector normal to the bounding sur- 
faces, £s = £o ^£r0s and (ps are the surface values of e and 
(/i, respectively. The last equation, just as for standard 
PB theory, expresses the electroneutrality of the system, 
as can be shown by the integral form of Gauss law. 

By solving the above set of equations, one can ob- 
tain the spatial profiles of the various degrees of freedom 
at thermodynamic equilibrium. For a general geometry, 
these equations can be solved only numerically. 



B. Bulk Behavior 

In the bulk, the system is homogeneous having a zero 
potential ip = and bulk values of n_|_ = n_ = rib and 



/)b- The EL equations, eqs [TTIIT3l reduce to: 

log(nba^) + Q;±0b - Ai± = 
log ( j + X(l - 20b) + a\a+ + a_)nb - = . 



(15) 



Eliminating the n± fields, we remain with a single bulk 
equilibrium equation 



log 



-t-x(l"20b)-f re-^("++° 



, 1 - f b , 
where F is defined as 

F = (a+ +a_)ei(''++'^-). 



= , 
(16) 

(17) 



Depending on the values of /i^ , /i± , a± , and X: the 
solutions of the bulk equation correspond either to a sin- 
gle phase of density (/)b or to a coexistence between two 
phases with different densities. Hereafter, we restrict our- 
selves to the single-phase region of the phase diagram, 
where the chemical potentials yU±, //^ follow directly from 
the form of the bulk free energy. 



C. Planar Geometry 

We exploit the symmetry of a planar system in order to 
derive analytically the pressure acting on the boundaries 
of the confined system. For a binary mixture confined to 
a slab delimited by two planar charged surfaces of infi- 
nite lateral extent (see Figure [1]), the general treatment 
introduced above can be simplified, and the free energy 
can be cast into a one-dimensional integral over the nor- 
mal z direction. For this special case, we show next that 
the pressure is proportional to the first integral of the EL 
equations. Using this expression, we also derive a first- 
order differential equation for the electrostatic potential 
that will greatly simplify the problem. 



1. Pressure in Planar Geometries 

We start from a general form of the free energy F which 
depends on N one-dimensional fields {'!/)i(z), ... , V'Ar(-z)} 
and their derivatives {ipiiz), ... , 7/1^(2)} 



F/A= / dzfi{^p,iz),^Pliz)};z) 



(18) 



When / docs not depend explicitly on the coordinate z, 
df/dz = 0, we obtain the following relation (see Ap- 
pendix): 



(19) 
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In our problem, / can be written as a sum of electrostatic 
and nonelectrostatic contributions 



N 



8tt 



i>'^ + ^ Qiniip + h ({rij}) , 



(20) 



where h is the grand potential of N different species with 
densities {ni, ... , titv} of a general form but without any 
electrostatic interactions. The charge of the ith species 
is denoted by qi, and e {{rii}) is the dielectric response 
as a function of the densities {ni, ... ,njv}. Substituting 
eq [20] into eq[T9l we obtain 



N 



const . 



(21) 



Finally, using the equilibrium equations for the densities 
{rii, ... ,n7v} 



Of 



1 de 



drii Stt drii ^ 



an,. 



(22) 



we end up with the following expression 



1 

8^ 



Eut 



Edh 



const . (23) 



For the special case of non-charged liquid mixtures, / re- 
duces to /ijwhile it follows from general thermodynamic 
identities [3^ that the last two terms in eq [221 are equal 
to the negative of the local pressure 



dh 



(24) 



total zz component of the stress tensor, which in ther- 
modynamic equilibrium has to be a constant and equal 
to -P, eq[25l 

Note that the above proof is valid for any form of the 
free energy / (as h had an arbitrary form), and accounts 
for electrostatic as well as non-electrostatic degrees of 
freedom in a completely general way. Applying this gen- 
eral result to our free energy, eqs[5][71 yields the following 
form of the total pressure: 



1 



P= -— (£o-2e,0)V 



/2 



Stt 



log(l 



+kBT ( a+n+(p + a^n-cf) ~~^^-^ \ . (26) 



This pressure P should be compared with the pressure 
of the standard PB theory: 

PPB = - -^£^'2 + fceT (n+ + 7i_ ) , (27) 
Stt 

and contains several additional terms. In fact, the dif- 
ference between the two is twofold: first, a polarizabil- 
ity term of the form Sr^'ip''^ /4,tt is included in the pres- 
sure since the dielectric constant is now spatially depen- 
dent. This term is equal to the product of the polar- 
ization —e^cfiip' and the electric field E = —i/j'- Second, 
the short-range interactions also change the form of the 
pressure: the solvent interactions contribute the term 
~X0^/fl'^! and the ion-solvent interactions add the two 
terms, a+n^cf) and a_n_(/). This last addition changes 
the pressure significantly when considering two similarly 
charged surfaces. We will discuss this point at length 
below. 



However, even in a charged liquid mixture the electro- 
static potential vanishes away from the boundaries so 
that P is also the bulk value of the pressure in a charged 
system. Together with eq [21] it follows that the first 
integral can be cast into the form: 



P=l- 
Stt 



Ede 



(25) 



Namely, the integration constant of eq [221 is simply the 
negative of the pressure, and is a constant throughout 
the system. We next consider separately the properties 
of the electrostatic and non-electrostatic terms in eq [221 
The first term is nothing but the negative of the zz 
component of the Maxwell electrostatic stress tensor, ap- 
propriately generalized to the case where the dielectric 
permeability is density dependent (33 |. The last two 
terms together, as already noticed, represent the local 
pressure of the system in the presence of charges. In 
the standard PB theory, these two terms are given by 
the van't Hoff form, while here they are given by an ap- 
propriate generalization, stemming from the free energy 
ansatz, eq[21 Combining all the terms in eg 1231 we get the 



2. First Integral of the EL Equations in Planar Geometries 

We now use the form of the first integral of the EL 
equations feq I26p to obtain an explicit first-order differ- 
ential equation for the electric field. The EL equations 
for the ion densities, eq[T2l give the following relations: 

n±iip, (/)) = „be=F«'^/'^'=^^-"±('^-'^'=) . (28) 

From the first integral, we now deduce 

d^y SirkBT ( 
dz ) (so - 2er4>) \ 

+a+{n+4) - (phnh) + a^{n^(l) - (?!)b"b) 
6b 1 



1 , 1 

-7 log 



1 - 



q^t) - nbB . (29) 



The difference between the pressure P at finite separation 
and its bulk value Pf, (infinite separation) is given by the 
rescaled osmotic pressure 11 = (P — Pb) /k-QTn\, and 



Ph 



2 + (a+ + a_)0b- 



1 



log(l-0b)- 



(30) 
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For a single plate (or, equivalently in the limit of two 
plates at infinite separation), 11 vanishes. Note that n± 
in eq[29]are functions of (f> and ip, and (j) is by itself a func- 
tion of ip and ?/;', given by eq [12] that is a transcendental 
algebraic equation for 0: 

/ d-0 \ ^ , , , , 

4(log^-log^-2,(^-<^,))=0.(31) 

The boundary conditions for each plate/boundary are 
given by three coupled algebraic equations for t/ig, ips and 
ipg. The first two equations are eqs [25] and [3TJ The 
third equation is given by the electroneutrality condition, 
eg I14[ that can be simply rewritten in the form 

Esips + 47recr = , (32) 

and Eg was defined after eg [HI 

III. RESULTS AND DISCUSSION 

The equilibrium equations eqs 1291 1311 and [35] derived 
above have no closed analytical solution. Hence, we solve 
them numerically to obtain spatial profiles for and n± . 
In addition, by considering the new terms as small per- 
turbations (compared to the regular PB theory), we show 
that an approximate analytical solution can be derived 
for the single plate case in the absence of salt. We then 
show numerical results for the pressure as a function of 
separation and for the pressure dependence on the exper- 
imentally controlled parameters a+, 0b, and n^,. Lastly, 



we compare our results to one available set of experiments 
on DNA in a binary solvent mixture. 



A. Density and Permeability Profiles 

We investigate the limit of small concentrations of the 
low dielectric component with no preferential solvation 
interactions {a± = 0) and for two values of the x param- 
eter. In Figure [21 we compare these numerical solutions 
for a single surface (for which the osmotic pressure van- 
ishes, n = 0, similar to infinite interplatc separation) to 
the ones of the regular PB equation with a homogeneous 
(average) dielectric constant, £av = £a — 4>h{£A — eb)- 

When all additional interactions are omitted {a± — 
and X = 0), the difference between the two models is 
negligible. While the deviation in cf) right at the sur- 
face reaches 10% of its bulk value 0b, it leads to only 
a 0.5% deviation for the dielectric constant e at the sur- 
face. The dependence of the other fields ip and n± on cj) is 
only due to changes in the dielectric constant. Therefore, 
in the limit of no short-range interactions, these fields 
hardly differ from the results of the regular PB model 
with homogeneous dielectric constant eav The correc- 
tion due to addition of solvent short-range interactions 
is also found to be small, even for larger solvent-solvent 
interaction, x = 1-5. This x value still describes a single 
bulk phase, as it is smaller than the critical value Xc = 2. 
We conclude that, in the absence of preferential solvation, 
a± = 0, the modified PB model has only a small added 
effect on the permeability, as can be seen in Figure [5] 



In Figure [3] we examine numerically the effect of pref- 
erential solvation on the solvent profile in the low concen- 
tration limit (0b = 0.09). Simply stated, when the ions 
prefer to be in the vicinity of the high permeability sol- 
vent molecules, we expect an increase in the exclusion of 
the low permeability solvent near the wall. Indeed, as can 
be seen in Figure [3] the exclusion of the low permeability 
solvent increases with a+ . We also find that the value of 
the co-ion-solvent short-range interaction parameter a_ 
does not significantly affect the form of the permeability 
profile and is set hereafter to zero. 

Moreover, when increasing a+ even further {a+ = 30 
in Figure[3]), (t>{z) and e{z) have a sharp change at about 
z = 2A. Namely, a layer rich in A species is formed near 



the wall with a thickness of a few angstroms, and at a 
certain distance from the wall decreases abruptly to 
a value close to the bulk value 0b. This phenomenon 
is clearly a consequence of the 7i+0 coupling and is in 
itself not an electrostatic effect. A gradient squared 
term (V0) in the free energy would have smoothed out 
this behavior and will be considered elsewhere. The 
steep variation observed for the nonhomogcneous mix- 
ture could be due to the fact that the boundary condi- 
tion demands a value that is much different from the 
one that follows from the bulk equation feg fTB)). so that 
the system tends to exhibit bulklike behavior as soon as 
possible. A detailed analysis of this phenomena will be 
presented in a separate study. 



B. Some Analytic Results ing certain limiting behaviors. We assume that there 

The model described above can be solved analytically 
by making some simplifying assumptions, and consider- 
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FIG. 2: Spatial profiles of (a) the solvent relative composition (j) and (b) the permeability e. The regular PB with homogeneous 
dielectric constant e = 77 (solid line) is compared with our modified PB for binary mixture with and without short-range 
interactions, x = (dashed line) and x = 1-5 (dotted line), respectively. Other parameters are: a = — 1/lOOA^^, rib = lO^^M, 
£a = 80, £b = 20 and 4>h = 0.05. In all the cases, no preferential solvation is considered, a± = 0. 
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FIG. 3: Spatial profiles of (a) the solvent composition (j> and (b) the dielectric constant e for various values of a+ as shown 
in the legend. The co-ion parameter is a_ — 0, the surface charge density is cr = — 1/100A~^, and the salt concentration is 
rib ~ 10~*M. The solvent dielectric constants are £a = 80 and £b ~ 20. The bulk value of (f> ^ (f'b = 0.09 . 



are no solvent short-range interactions [x ~ 0) and that 
the preferential solvation interaction is weak compared 
to (a+ ^ 1). We also assume that the contribution 
due to the permeabihty inhomogeneities is neghgible, and 
take e throughout the system to be the weighted average 
of the two bulk relative compositions, £av Moreover, 
we take the limit of zero salt, as if only counterions are 
present to keep the system neutral. Lastly, we assume 
that one of the solvent concentrations is much smaller 
than the other. 

With these assumptions, we can practically isolate the 
effect of preferential solvation and obtain analytical pro- 



files. The PB equation in this limit assumes the form 
d^-iA _ 47ren(z) 

where the ion density n{z) is a function of both the po- 
tential -0 and the solvent relative composition (j) 

n{z) = nbe-ei'/kuT-a{4,-M^ (34) 

The prefactor rib is determined by satisfying the elec- 
troneutrality condition, and the subscript ± in n± is 
omitted in this counterion only case. 
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The composition as a function of n is 



9be 



(35) 



Here, we made use of the assumption that the preferen- 
tial solvation interaction is small, a{a^n) <^ 1. Substi- 
tuting it back into eg I34[ we obtain for the further limit, 
a(?i -Cl , 



ribC 



1 - Q;2a3(/ibnbe-'^''^/'=BT- 



(36) 



For finite values of eip/kBT and in the limit of 
a^cfiha^rih ^ 1, we obtain to lowest order 



n ~ TT-bC 



-eii/kgT 



l + . (37) 



Using this equation in the PB equation, we get a second 
order differential equation for 



(Pip iiren^- 



(38) 

For a single plate, the boundary condition to the equation 
above is given by 



A, 



GC 



(39) 



where = e^/eavfcsT', and Aqc is the well known Gouy- 
Chapman (GC) length 

Solving the equation above, we obtain the following 
results: 



^(z) 
n(z) = 



fcRT 



log [(z + ASc)'-A^] 



+ Age)' + \' 
27tIb [{z + ASc)2 - \l 



^0 , (40) 
(41) 



where Ag = a? 4)h / '^t^Ib is a typical length associated 
with the a parameter. The second length scale is the 
modified Gouy-Chapman (GC) length Aqq obtained by 
satisfying the boundary condition: 



A; 



A, 



GC 



GC 



(l + Vi + 4(A./Agc)2) . (42) 



For As < Aqc one obtains Ag^. ~ Aqc [l + (As/Aqc)^]- 
Thus, the effect of preferential solvation enhances the ion 
density in the proximity of the surface and results in a 
faster decay of the density profile. 



C. Pressure vs. Separation Curves 

Due to symmetry, the pressure between two identically 
charged plates is most conveniently calculated at the mid- 
plane. When considering the change in pressure due to 
the permeability inhomogeneity, one can separate the di- 
rect and indirect corrections. The direct one is due to 
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FIG. 4: Dependence of pressure on separation D for various 
ion-solvent interaction strengths a+ as shown in the legend. 
Other parameters are a = — 1/100A~'^ , rib = lO^^M, 
£A = 80, £B = 4, and 0b = 0.09 . 



the change of midplanc composition 4){z), while the in- 
direct one is related to changes of the midplanc ion den- 
sity. In the absence of preferential solvation {a± = 0), (/> 
at the midplanc depends only on the local value of the 
electrostatic field (see eq [3T] with a± ~ 0). However, 
in a symmetric twoplate system, the electrostatic field 
vanishes at the midplane and (j) there equals to its bulk 
value, ultimately contributing no correction to the pres- 
sure. Moreover, since e{z) turns out to be nearly homo- 
geneous (see Figure [2]) , the indirect correction is minute 
as well. Thus, in the absence of preferential solvation the 
combined effect of a binary mixture is negligible. 

When adding the preferential solvation term character- 
ized by the parameter a± (the term coupling between n± 
and 0), (j) becomes dependent on the nonzero midplane 
potential. As a consequence, the midplane (p value differs 
from (/)b , and results in two direct corrections to the pres- 
sure. The first comes from the osmotic pressure of the 
solvent k^T {(f) — (ph) / a^) , while the second originates 
from the (p ■ n± coupling term. Thus, even if the indi- 
rect contribution to the profiles is negligible, the direct 
correction alone can substantially change the pressure. 

The effect of preferential solvation on the pressure be- 
tween two identically charged surfaces is now examined 
for various values of the q;+ parameter. As the coupling 
term in eq [251 contributes directly to the pressure, the 
effect of changing a+ is rather pronounced in the pres- 
sure versus separation curves (see Figure U). The figure 
clearly shows that the pressure decreases when a_|_ in- 
creases. Yet, because the effect is short ranged, the major 
differences are observed at small separations {D < 2nm). 
From the analysis of the profiles with respect to the inter- 
action strength (Figure 2]) , we conclude that the change 
in pressure is substantial only up to distances of a few 
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FIG. 5: Pressure as a function of for various interplate 
separations: D — 0.5,1, and 2.5 nm. Other parameters are 
a = -I/IOOA"^ , rib = 10"^M, £a = 80, eb = 4, and 
0b = 0.09 . 



FIG. 6: Pressure versus bulk solvent composition 0b for var- 
ious ion-solvent interaction strength, a+. Other parameters 
are a = -1/lOOA"^ , Tib = lO'^^M, ea = 80 , and £b = 4 . 
The separation is fixed at D = Inm. 



nanometers. This means that for large separations the 
midplane values of the fields ijj, n and (j) within our model 
will be similar to the regular PB theory predictions. How- 
ever, for small separations of the order of 1-2 nanome- 
ters, both the profiles and the pressure are affected by 
the preferential solvation. 

D. The Effect of the Solution Parameters on the 
Pressure 

We proceed by examining the influence of the param- 
eters ck-i-, 0b and nb on the pressure. The a+ parameter 
can be modified experimentally by using different sol- 
vents, while 0b and nb can be easily controlled in the 
experiment by changing composition. 

In Figure [5l we present the dependence of the pressure 
n(Q!-|-; D) on the interaction strength a+ for a fixed sep- 
aration. As expected, for small values of a+ (< 2), the 
pressure depends only weakly on the interaction strength, 
whereas for larger values the pressure falls with a+ . This 
implies that there is a value of a+ where its direct con- 
tribution to the pressure becomes larger than all other 
contributions (electrostatic and entropic). Moreover, the 
slope of Il{a^; D) depends on the separation, as can be 
clearly seen by comparing the D = 0.5 nm and D = 1 nm 
results. For smaller separations {D = 0.5 nm), the pref- 
erential solvation effect is stronger in accordance with 
the results shown in the previous sections (Figures [2]-[4]) , 
where the effect has a range of a few nanometers. For 
D = 2.5 nm, the pressure changes very slowly with 
and the preferential solvation is small even for large 
(15 <a+< 20). 

Next, we investigate how increasing 0b changes the 
pressure. In Figure [51 we plot the pressure versus 0b 



for a fixed separation D = I nm for three values of the 
interaction strength Increasing the low permeabil- 
ity solvent concentration decreases the pressure through 
decrease of the permeability and by increasing the pref- 
erential solvation. The results in Figure [6] suggest that 
even for a+ (no preferential solvation) the pressure 
decreases with 0b, implying that the dielectrophoretic 
mechanism contributes a nearly linear dependence on e^- 
When increasing the n(0b) slope is steeper due to 
higher bulk pressure of the low permeability solvent, con- 
tributing directly to the pressure through the preferential 
solvation term feq [26)) . 

Since there is a linear 0b term in the pressure 
(~ —kBT(j>yj/a^), one can deduce from the results above 
that the main contribution to the pressure comes simply 
from a higher reference concentration 0b that reduces 
the pressure. Namely, for small separations, the solvent 
4> and the ion n± at the midplane have only a weak de- 
pendence on 0b, similar to regular PB theory where the 
pressure has no dependence on the bulk salt density nb 
at small separations [l3| . 

Finally, we investigate the influence of the salt concen- 
tration on the pressure at fixed separation {D = Inm). 
The results are presented in Figure 7, where we plot the 
pressure versus the salt concentration for fixed separa- 
tion and for different values of the interaction strength 
Q!-i_. For low concentration (nb < 10~^M), the pressure 
has no dependence on nb. It is known from the regular 
PB theory [l3| that, at small separations, when the De- 
bye length is much larger than the separation (Ad 3> D), 
the pressure only weakly depends on the salt concentra- 
tion. In this sense, the modified PB theory presented 
here is similar to the regular PB theory. The effect of 
preferential solvation is just a constant addition to the 
pressure, which keeps the same dependence of the pres- 
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FIG. 7: Pressure versus bulk salt concentration Wb for various 
ion-solvent interaction strength, a±. Other parameters are 
a = -1/100A~^ EA = 80, Eb = 4, and = 0.09. The 
separation is fixed at _D = Inm. Lines are guides for the eye. 



sure on the salt concentration. When going to larger 
concentrations, where the Debye length becomes com- 
parable to or smaller than the separation, the pressure 
decays exponentially as can be expected also from stan- 
dard PB theory. 



E. Comparison of Model with Experiment 

So far, we have shown that we are able to account 
for some of the forces that lead to cosolvent preferen- 
tial exclusion from charged interacting macromolccular 
surfaces. In particular, our model accounts not only for 
locally varying dielectric profiles that follow the solvent 
mixture composition (through the variables Ea and £b), 
but also for preferential ion solvation (through the vari- 



ables Q;±), that in turn depends on local solvent com- 
position as well. Using these two sets of parameters, 
it is possible to propose a physical mechanism for so- 
lute (or solvent) exclusion from interacting surfaces. The 
different dielectric constants of the two solvent compo- 
nents cause depletion of one of the components from the 
charged surface. This variation in solvent composition 
can in turn affect the concentration of ions between the 
two interacting surfaces. The combined effect can lower 
the disjoining pressure between equally charged surfaces 
by varying local solution concentrations. 

While it is hard to unambiguously prove the origins 
of the molecular interactions that lead to the different 
solvation properties of ions and cosolvent in the vicin- 
ity of complex macromolccules, we show that by varying 
the model parameters we can account for the observed 
trends in the experimental studies of Rau and Stanley 
@. In these experiments, osmotic pressure is applied 
to a condensed phase of DNA strands in aqueous solu- 
tion by adding a neutral polymer, poly(cthylenc glycol) 
(PEG) that is completely excluded from the DNA phase. 
The DNA-DNA spacings in solution, D, are measured 
using small angle X-ray scattering. In addition, salt and 
different alcohols are added at different concentrations. 
Figure [5^ shows the experimentally derived 11(1?) (equa- 
tion of state) for DNA solutions containing either 0.02 
or 1.2 M NaBr salts, as they appear in rcf [6|. Solu- 
tions to which 2-methyl-2,4-pentanediol (MPD) alcohol 
was added at concentrations of 0.5 or 1 M are compared 
with solution with no MPD. As the figure shows, for both 
salt concentrations, the added salt lowers the DNA-DNA 
spacing for a particular applied osmotic pressure. How- 
ever, the figure also clearly shows that the reduction in 11 
is more significant at the higher salt concentration. This 
would imply that salts (or more generally electrostatic 
forces) are involved in determining the effect of MPD on 
n, suggesting an important role reserved for the dielectric 
properties that can be linked to the distribution of ions 
and cosolvent partitioning in between the DNA strands. 



Even though our model considers interactions between 
flat surfaces, rather than the cylindrical ones expected 
for DNA strands, we show that it is possible to get sim- 
ilar trends to those found in experiment using the two 
sets of the e and a parameters. Figure [5)d shows our re- 
sults for two charged plates with charge density similar to 
that of DNA (tr = — 1/100A~^), and dielectric constants 
of tA = 80 (representing water) and = 25 (close to 
the value of pure MPD). In our model, we use only the 
relative volume fraction of the two species, water and 
MPD, without accounting for their different molar weight 
(water molar weight is 18 and that of MPD is 118). In 
the experiment [0], the two MPD solutions have concen- 
trations of 0.5 and 1.0 M corresponding, respectively, to 



4>\, = 0.055 and 0.11. These values are in good agree- 
ment with the values chosen in our model (Figure ^p) 
to give a good fit to the experimental data: (p^ = 0.126 
and 0.252. Note that the value (j)h = 0.0018 is chosen 
for convenience to fit the zero MPD concentration case. 
The value of a+ could, in principle, be evaluated from 
the solvation free energy of NaBr salts in binary water- 
MPD solutions of different contents. Because such data 
is lacking, we treat a+ as a fitting parameter and use 
a+ = 5 (in units of kuT), allowing us to fit closely the 
experimental data. We note that this a+ value is close 
to the experimental transfer free energy of sodium from 
water to cthanol ~ bk^T as was quoted in rcf [s^]. The 
value of a_ for bromide, which was taken in our work as 
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FIG. 8; Pressure IT as a function of the separation D. (a) Experimental data. Circles (squares) represent the 0.02 M (1.2 M) 
added NaBr salt results adopted from ref Empty circles and squares represent the experiment with no MPD. Gray circles 
and squares represent the experiment with 0.5 M MPD. Black circles and squares represent the experiment with IM MPD. 
The dashed lines are guides to the eye. (b) Numerically calculated lines from the model with parameters taken to match the 
experiment. The ion-solvent interaction strength is the same for all lines a+ = 5 , and it was treated as a fitting parameter. 
The surface charge is taken as a = — 1/100A~^ to fit the DNA values. The dielectric constants are ea ~ 80 (water) and 
£b = 25 (MPD). Black lines are for rib — 0.02 M and gray lines are for rib = 1.2 M. Solid, dashed, and dotted lines stand for 
0b = 0.0018 , 0b = 0.126 , and 0b = 0.252 , respectively. 



zero, seems to be generally lower than that of but is 
less clearly resolved [Slj . 

We find that for this set of parameters the numerical 
calculation grosso modo follows the experimental trends: 
the spacing between curves grows with equal additions of 
alcohol to the solution, but for the higher salt concentra- 
tion the change in 11 is larger. These results underscore 
two important general conclusions. First, it would be im- 
possible to explain the difference in the n(D) behavior 
for high salt versus low salt without discussing ions and 
the role of the dielectric medium. Our model introduces 
these species in a self-consistent manner through the PB- 
like theory. Second, the comparison demonstrates the 
important role reserved for the preferential ion-solvation 
interactions in the different components of the binary so- 
lution. Specifically, it would be impossible to explain 
the shifts in distances at a given applied osmotic stress 
without invoking a nonzero a+, that in turn shifts the 
exclusion of ions due to changes in solution composition 
in the DNA phase. 



IV. CONCLUSIONS 

The model presented in this work is a modification of 
the regular PB theory and accounts for the effects of a di- 
electric medium composed of two solvents. There are two 
important features that modify the standard PB theory. 
Due to different dielectric constants of the two solvents, 
the permeability is no longer homogeneous. This effect 



is accounted for by coupling the electrostatic field to the 
solvent local composition. The second modification is the 
addition of a preferential solvation term which enables 
the ions to drag a favorable solvent and locally deplete 
the other. This is modeled by a coupling term between 
the ion density and the solvent local composition. Simi- 
lar models can be used for a number of problems such as 
the properties of interfaces [l^, [l^ , critical behavior of 
binary mixture in the presence of ions [l^, [l3| , and the 
forces between more elaborated charged macromolecules, 
which is the main focus of this work. 

For two identically charged planes, we find that, in 
the absence of preferential solvation, the density pro- 
files and the pressure undergo only small modifications. 
This is demonstrated numerically and supported by an 
argument that relies on the symmetry of the two-plate 
system. However, by adding the preferential solvation 
term, we are able to observe a considerable correction 
to the pressure at small separations. The coupling be- 
tween ion density and solvent local composition appre- 
ciably changes the midplane concentration values, and as 
a consequence the pressure is reduced. We also investi- 
gated the dependence of the pressure on experimentally 
controlled parameters such as salt concentration, bulk 
solvent composition, and preferential solvation strength. 
The pressure depends on the preferential solvation but 
changes substantially only at small separations (1-2 nm). 
The threshold of preferential solvation energy which is 
required to change the pressure significantly is on the or- 
der of a few k^T. It is found that the pressure depends 



12 



nearly linearly on the bulk relative composition, imply- 
ing that the density profiles at small separations have no 
dependence on 0b- Finally, the effect of added salt to 
the solution seems to change the pressure in the same 
manner as for regular PB theory. 

We also used our model to put fourth an explanation 
for the experimentally measured pressures in a condensed 
phase of DNA. The comparison shows that the experi- 
mental trend is bourne out by our model results. This 
suggests that the main mechanism causing the deple- 
tion of one solvent away from the charged macromolecule 
is very plausibly the preferential solvation of the ions. 
Thus, beyond simple electrostatic screening, salt ions 
may play an additional and important role in the be- 
havior of charged macromolecules immersed in solution. 

Further applications and refinements of the model 
could be considered. For example, the model can be 
used to analyze the effect of strong preferential solva- 
tion on the critical behavior. As shown in Figure [H the 
solvent relative composition profile becomes discontin- 
uous at strong preferential solvation. We believe that 
this can be explained in the framework of a Ginzburg- 
Landau theory that would account for the phase transi- 
tion. Moreover, in addition to the simple planar geom- 
etry treated here, other geometries such as a cylindrical 
one can obtain a more direct quantitative comparison 
with DNA experiments. Finally, the limit of ionic dilute 
solutions can be generalized to the concentrated limit, 
including the full entropy of mixing, as was considered in 
ref [si. 
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APPENDIX A: DERIVATION OF PRESSURE IN 
ONE DIMENSIONAL SYSTEM 

In the following we present the derivation of eq. [TH 
We start from a free energy F, eq 18 which depends on 
N coupled one-dimensional fields {^i{z), ... ,'iIjn{z)} and 
their derivatives {ip[{z), ... ,ip']y{z)}: 



F/A= / f{{Mz),i'liz)};z)dz, 



where i = 1, ... , A^. There are N EL equations 
SF 



Hi 



= 



dlpi dz dtp[ 



The total derivative of /(z) is 



Hr I Pl,l„ ^» dtp- -^r 



dz ^ V dth. 

i ^ 

Moreover, we can write 



dz 



dz\diP'J'J dz 



(Al) 



(A2) 



(A3) 



(A4) 
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Substituting this back into eq |A3l and using the EL equa- 
tions we find 



dz 



E 



df 

dip: 



dj_ 

dz 



(A5) 



When /(z) does not depend explicitly on the coordinate 
z, ^ = 0, the last term vanishes and we end up with a 
first order differential relation: 



/-E 



df 



dtp. 



7^1 



const . 



(A6) 
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